Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 7 October 2009 



(MN MfeX style file v2.2) 



The Imperial IRAS-FSC Redshift Catalogue: luminosity 
functions, evolution and galaxy bias 

Lingyu Wang 1 ' 2 *, Michael Rowan- Robinson 1 

1 Astrophysics Group, Blackett Laboratory, Imperial College of Science Technology and Medicine, London SW7 2AZ 

2 Astronomy Centre, Department of Physics and Astronomy, University of Sussex, Falmer, Brighton, BN1 9QH 



Accepted . Received ; in original form 



ABSTRACT 

We present the luminosity function and selection function of 60 fim galaxies selected 
from the Imperial IRAS-FSC Redshift Catalogue (IIFSCz). Three methods, including 
the 1/Vmax and the parametric and non-parametric maximum likelihood estimator, 
are used and results agree well with each other. A density evolution oc (1 + z) 3 - 4±0 9 or 
a luminosity evolution oc exp(1.7 tj^/r) where ti, is the look-back time is detected in 
the full sample in the redshift range [0.02, 0.1], consistent with previous analyses. Of 
the four infrared subpopulations, cirrus-type galaxies and M82-type starbursts show 
similar evolutionary trends, galaxies with significant AGN contributions show stronger 
positive evolution and Arp 220-type starbursts exhibit strong negative evolution. The 
dominant subpopulation changes from cirrus-type galaxies to M82-type starbursts at 
logioCZW^o) ~ 10.3. 

In the second half of the paper, we derive the projected two-point spatial correla- 
tion function for galaxies of different infrared template type. The mean relative bias 
between cirrus-type galaxies and M82-type starbursts, which correspond to quiescent 
galaxies with optically thin interstellar dust and actively star-forming galaxies respec- 
tively, is calculated to be 6 c irrus/^M82 = 1-25 ± 0.07. The relation between current 
star formation rate (SFR) in star-forming galaxies and environment is investigated 
by looking at the the dependence of clustering on infrared luminosity. We found that 
M82-type actively star-forming galaxies show stronger clustering as infrared luminos- 
ity / SFR increases. The correlation between clustering strength and SFR in the local 
Universe seems to echo the basic trend seen in star-forming galaxies in the Great 
Observatories Origins Deep Survey (GOODS) fields (at z ~ 1). 

Key words: galaxies: evolution - galaxies: statistics - galaxies: luminosity function, 
mass function - infrared: galaxies - large-scale structure of Universe. 



1 INTRODUCTION 

The study of galaxy distribution on large scales (~ 10 ft -1 
Mpc) where density fluctuations are still in the linear regime 
(Sp/p <C 1) has been an important tool in constraining cos- 
mological models. However, the unknown relationship be- 
tween the distribution of dark matter and that of luminous 
objects, usually referred to as galaxy bias, has complicated 
the extraction of cosmological parameters, especially when 
we compare results from surveys of different galaxy subpop- 
ulations or of varying depths. At small scales (~ 1 h^ 1 Mpc) 
where nonlinear effects become important, detailed studies 
of galaxy bias can improve our understanding of galaxy for- 
mation and evolution which is still unclear due to the con- 
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volution of gravitational evolution and complex astrophys- 
ical processes (e.g. gas cooling, fragmentation, star forma- 
tion, feedback) . This paper attempts to measure the relative 
galaxy bias over the intermediate scale (1 ~ 10) Mpc, i.e. the 
weakly nonlinear scale, as we are constrained by low galaxy 
density at small scales. 

Several models of galaxy bias have been developed over 
the past 20 years or so. The simplest one is to assume galax- 
ies trace the underlying mass in a linear, deterministic way, 
5 g — bS m where 5 S and 8 m is the galaxy and matter density 
perturbation (or contrast) field and b is the linear bias pa- 
rameter. In the peaks bias model (Kaiser 1984; Bardeen et 
al. 1986), the relationship between the correlation function 
of galaxies and that of the underlying dark matter satisfies 
£g( r ) ~ fr 2 £m(y), based on statistics of the density peaks 
in a Gaussian random field. On large scales, it seems rea- 
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sonable to assume a scale independent linear galaxy bias if 
galaxies form in virialized dark matter halos. On interme- 
diate and small scales, the situation is more complicated. 
In more realistic galaxy bias models, b is thought to be 
time-varying, scale-dependent, non-local, or even stochas- 
tic to some extent. N-body simulations of dark matter and 
semi-analytic models of galaxy formation predict a scale- 
dependent relation between galaxies and dark matter due to 
the mass- dependence of galaxy formation efficiency as well 
as a luminosity-dependent bias for galaxies with luminosities 
above L* (Benson et al. 2000; Benson et al. 2001). Attempts 
to understand galaxy bias have also been made in the halo 
occupation distribution (HOD) framework where the rela- 
tionship between the galaxy distribution and the underlying 
dark matter is encoded in the conditional probability func- 
tion P(N\M) which gives the number of galaxies of a given 
type in a halo of virial mass M (Jing et al. 1998; Ma & Fry 
2000; Peacock & Smith 2000; Cooray & Sheth 2002; Zehavi 
et al. 2005; Tinker et al. 2005; Zheng, Coil & Zehavi 2007; 
Simon et al. 2008). 

Observationally, several manifestations of galaxy bias 
are known. From dipole analysis, different values of the 
ft = Q^ 6 /6 parameter have been reported using different 
types of mass tracer (Rowan-Robinson et al. 2000; Ko- 
cevski, Mullis & Ebeling 2004; Erdogdu et al. 2006 and refer- 
ences therein) . The well-known morphology- density relation 
(Dressier 1980; Goto et al. 2003) states that late-type galax- 
ies (disk-dominated) prefer low density regions and early- 
type (bulge-dominated) galaxies dominate in dense cluster 
cores. There is also a correlation between colour and mor- 
phology (the majority of blue galaxies are late-type galaxies, 
while red galaxies are mostly early- type galaxies), as well as 
a correlation between colour and luminosity (red galaxies are 
systematically more luminous than blue galaxies). There- 
fore, galaxy clustering is expected to be dependent on phys- 
ical properties such as morphology, colour, surface bright- 
ness and luminosity. Using a total of over 100, 000 galaxies 
from the 2dF Galaxy Redshift Survey (2dFGRS), Norberg et 
al. (2001) found a scale-independent luminosity dependence 
of galaxy bias of the form b/b* = 0.85 + 0.15L/L* which 
means the luminosity-dependence is conspicuous above V . 
It agrees with the conclusions of Benoist et al. (1996) based 
upon 3,600 galaxies selected from the Southern Sky Red- 
shift Survey 2 at the faint end, but it exhibits a less steep 
rise at the bright end. Zehavi et al. (2002, 2005) studied 
galaxy clustering in the SDSS, which not only confirmed 
the conclusions of Norberg et al. (2001) but also probed 
the relative bias behaviour to fainter galaxies. Madgwick et 
al. (2003) divided their sample from the 2dFGRS into rel- 
atively passive and active star-forming galaxy subsamples. 
A scale-dependent relative bias was found in comparing the 
clustering strength of the two subsamples. On scales smaller 
than 8 Mpc, the relative bias between the two spec- 
tral classes is & pa ssivc/&activc = 1-45 ± 0.14. On scales larger 
than 10 /i -1 Mpc, the relative bias reduces to unity. Zehavi 
et al. (2002; 2005) used a luminosity- dependent colour-cut 
to divide their volume-limited sample from the SDSS. The 
colour-dependence of the correlation function is found to 
be very similar to that found in Madgwick et al. (2003). 
In addition, the luminosity- dependence for red galaxies is 
such that both faint and luminous red galaxies cluster more 
strongly than moderately bright red galaxies around L* 



while a luminosity-dependence for blue galaxies is not no- 
ticeable (Wang et al. 2007; Swanson et al. 2008; Cresswell 
& Percival 2009). 



The optical colour of a galaxy is connected with a com- 
plicated combination of star formation history and current 
star formation rate (SFR). In contrast, it is well known that 
the far-infrared emission is a good SFR indicator due to 
the re-processing of starlight by dust. So in principle one 
can measure clustering strength of infrared galaxies to learn 
about the relation between environment and SFR. So far, 
little is known about the clustering of infrared galaxies and 
its dependence on various physical parameters. The largest 
sample of local infrared galaxies is provided by IRAS. The 
correlation length of IRAS galaxies (ro ~ 4 ft -1 Mpc) is 
similar to that of local star-forming/blue galaxies found in 
the SDSS and 2dFGRS (Saunders et al. 1992; Fisher et al. 
1994). Clustering measurements of infrared galaxies did not 
show significant dependence on infrared luminosity (Beis- 
bart & Kerscher 2000; Szapudi et al. 2000; Hawkins et al. 
2001). This was thought to be a result of the combined ef- 
fect of the dark matter halo mass and SFR in determining 
the far-infrared luminosity. From the SFR - density relation, 
we know that star-forming galaxies preferentially reside in 
low density environments (Gomez et al. 2003). On the other 
hand, the gas content in a galaxy is largely dependent on the 
mass of the host dark matter halo. The dependence of clus- 
tering on infrared galaxy types is not entirely clear either. 
With a few thousand galaxies in the QDOT IRAS galaxy 
redshift survey (Lawrence et al. 1999), Mann, Saunders & 
Taylor (1996) did not detect any significant variation in the 
clustering of warm and cool IRAS galaxies separated by the 
far-infrared emission temperature (above or below 36 K re- 
spectively), assuming the far-infrared spectrum can be fitted 
by S v oc vB v (T\ Hawkins et al. (2001), using the much big- 
ger IRAS PSCz sample (Saunders et al. 2000), found that 
the relative bias between cool and warm infrared galaxies is 
around 1.5 on scales between 1 and 10 hT 1 Mpc. 

Our aim is to investigate the relation between star for- 
mation activity in star-forming galaxies and environment 
using infrared samples in the general field in the local Uni- 
verse. We study the clustering dependence on infrared lumi- 
nosity and infrared template type, taking advantage of a new 
redshift catalogue of over 60,000 galaxies selected at 60 /im 
from the IRAS Faint Source Catalog (FSC). The remain- 
der of this paper is organised as follows. Section [2] briefly 
describes the galaxy sample used. Details of the new cata- 
logue can be found in Wang & Rowan- Robinson (2009). In 
Section [3] we derive the luminosity function and the selec- 
tion function using various methods which will then be used 
to construct random catalogues in measuring the clustering 
signal with flux-limited samples. In Section [4] clustering de- 
pendence on infrared template type and infrared luminosity 
is investigated using the projected two-point spatial corre- 
lation function. Finally, discussions and conclusions are pre- 
sented in Section [5] Throughout the paper, unless otherwise 
stated, we assume Ho = 72 km s _1 Mpc -1 and a flat uni- 
verse with SIm = 0.3 and £Ia = 0.7. Luminosities are quoted 
in units of Lq. 
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2 THE GALAXY SAMPLE 

The Imperial IRAS-FSC Redshift Catalogue (IIFSCz; Wang 
& Rowan- Robinson 2008) is a recently constructed cata- 
logue containing 60,303 galaxies selected at 60 /im from the 
IRAS FSC. The IIFSCz covers most of the sky at \b\ > 20° 
(~ 61% of the entire sky). 55% of the sources have spectro- 
scopic redshifts obtained from past redshift surveys and a 
further 20% have optical, near- infrared and/or radio iden- 
tifications for which photometric redshifts were estimated 
through the training set or template-fitting method. We also 
determine the best-fit infrared template for each galaxy with 
either spectroscopic or photometric redshift. Fig. [T] shows 
the main four infrared template types, infrared cirrus (dust 
emission due to the interstellar radiation field), M82 star- 
burst (interaction-induced starburst), AGN dust torus emis- 
sion and Arp 220 starburst (merger- induced starburst). 

The intrinsic 90% completeness limit of the IRAS 
FSC is at around F^o = 0.36 Jy. In order to reduce 
incompleteness-induced problems, we have selected a sam- 
ple of around 22,500 galaxies with Fao 0.36 Jy. At this 
flux limit, around 75% of our sample have spectroscopic red- 
shifts, 15% have photometric redshifts and 10% do not have 
any redshift estimates. Almost all galaxies in areas covered 
by the SDSS DR6 survey have redshift, either spectroscopic 
or photometric. Thus, the redshift completeness is a spa- 
tially varying function and this is taken into account in this 
paper. 

We convert heliocentric redshift to redshift in the Local 
Group frame. Redshifts are used as distance indicators with- 
out corrections for peculiar motions. In consequence, we do 
not include galaxies below redshift z — 0.003 for which z is 
a poor estimate of the distance, in our sample. In Fig. the 
redshift distribution of galaxies in our sample is plotted in 
comparison with galaxies in the IRAS PSCz. Fig. [3] shows 
the normalised redshift distribution for each infrared galaxy 
type in our sample. 

Cirrus-type galaxies are mainly found in the local Uni- 
verse (z < 0.06), while the other three types are more com- 
mon at relatively high redshift. So, cirrus-type galaxies are 
almost absent at the high luminosity end. This transition 
in the predominant subpopulation in the mixture across the 
redshift (or luminosity) range is important in understanding 
the clustering dependence on luminosity for the full sample. 
We will return to this point in Section 4. 
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Figure 1. The four infrared templates used to fit IRAS fluxes 
from 12 to 100 fim for galaxies in the IIFSCz. 
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3 THE LUMINOSITY AND SELECTION 
FUNCTION 

The selection function S(z) (hereafter SF) of a flux-limited 
survey is defined as the mean comoving number density of 
galaxies that can be observed at redshift z. It is the inte- 
gration of the luminosity function $(-£/) (hereafter LF) at 
redshift z which specifies the luminosity distribution of the 
galaxy population observed in the survey, 



S(z) 



$ z (L)dL, 



(1) 



where L m i n ( z ) is the minimum luminosity a galaxy at redshift 
z can have in order to remain above the flux limit. The 
inherent assumption is that the LF does not depend on the 



Figure 2. The redshift distribution of the IRAS PSCz galaxies 
(thin solid line) and IIFSCz galaxies with either spectroscopic or 
photometric redshifts (thick solid line) or IIFSCz galaxies with 
spectroscopic redshifts only (dotted line). The dashed curves are 
obtained from the selection functions of the PSCz (thin dashed 
curve) and IIFSCz (thick dashed curve). Galaxies with redshifts 
z < 0.003 are excluded. 



local density field. However, it has been known for some 
time that the LF is dependent on galaxy type which is in 
turn related to the environment. Therefore, the assumption 
that the bivariate distribution function of luminosity and 
position is separable is not strictly valid. 
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Figure 3. The redshift distribution of each infrared population 
in the IIFSCz. The peaks are normalised to 1. 



3.1 Methods 

The LF is of fundamental importance in constraining mod- 
els of galaxy formation and evolution. Various techniques 
have been proposed to measure it, such as the maximum- 
volume estimator 1/Vmax (Schmidt 1968) and the paramet- 
ric (Sandage, Tammann & Yahil 1979, hereafter STY) and 
non-parametric maximum likelihood estimator (Efstathiou, 
Ellis & Peterson 1988; Springel & White 1998, hereafter 
SW98). The 1/V ma x method assumes a uniform number den- 
sity throughout the observed volume and therefore is vul- 
nerable to density inhomogeneities present in a survey. The 
maximum likelihood estimators have the advantage of being 
insensitive to any local clustering effect, but lose information 
on the overall normalisation. 

The STY method requires a prior knowledge of the 
appropriate functional form of the LF. The probability 
of a source at redshift Zi having luminosity in the range 
[Li , Li + dL] is 



p(Li\zi) 



n(zi,Lj) _ $ Zi (Li) 

/~. n(zi,L>)dL> ~ S{z % ) 



(2) 



In Bayesian statistics, the posterior probability of the pa- 
rameters characterising the LF given the data is propor- 
tional to the likelihood function (i.e. the probability of the 
data given the parameters) assuming a uniform prior prob- 
ability of the parameters. Therefore, the best-fit parameters 
in the LF are the ones which maximise the likelihood given 
by 



C = Y[p{Li\zi) 



Figure 4. The x-axis is the 60 fim luminosity with a type- 
dependent K-correction derived from the best-fit infrared tem- 
plate. The y-axis is the 60 fim luminosity with a K-correction 
derived from a power-law SED S„ cx u~ 2 . The dotted line corre- 
sponds the perfect agreement between the two. The blue dots are 
AGN dust tori and the black dots are galaxies of any other type 
in the redshift range 0.003 ^ z < 4.0. 



where C = 2.6 x 10~ 2 h 3 Mpc" 3 , a = 1.09, a = 0.724, 
L* = 10 8 ' 47 h~ 2 L Q . Note that ip{L) is defined as the LF 
per decade in luminosity. The normalisation C is derived 
from the expected number density of sources with Fm J? 0.6 
Jy per unit solid angle. 

The parametric maximum likelihood technique does not 
have an intrinsic measure of goodness-of-fit. Therefore, a 
non-parametric approach is needed to justify the functional 
form of the LF or the closely related SF. SW98 proposed a 
non-parametric maximum likelihood estimator which adopts 
stepwise power laws to characterise the SF, 



S(z) =S k [ — 



(5) 



where k is the redshift bin number. For redshift in the k bin 
(xh-i < z ^ xt), rrik is the logarithmic slope of the power 
law. 

If we assume a pure density evolution of the form 

n(z) = n g(z) = n (l + z) P , (6) 

the best estimates of mt can be derived from the following 
equation, 

d(ln£) Sk, ai 



(3) 



dm k 
where 



+ T k = 0. 



(J) 



Saunders et al. (1990; hereafter S90), using the STY 
method, determined the present-epoch 60 /im LF from 2818 
IRAS galaxies which is a Gaussian combined with a power 
law, 



ip(L)=Cl — 



exp 



1 1 2 



(4) 



2. 



^— - + 5k, ai In x 



The variance of m k has the simple form 



In- 



.(8) 



V[m k ] = < y^<?fc,q 



rn k - z t 



g'izT) 



(9) 
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Note that z™ is the maximal redshift a galaxy can have such 
that L m i n (zY i ) = Li. Without including luminosity evolu- 
tion, z™ can be derived by solving the following nonlinear 
equation, 



Fj l + zi 
F min 1 + z? 



(10) 



The comoving distance in a ACDM cosmology can be cal- 
culated from 



dz_ 

Ho J y/siuO- + z') 3 + n A ' 

Here $ denotes the K-correction 
/ R{y)Uv n )d V 



/ R{v')f v {v')dv' 



(11) 



(12) 



where R(u) is the responsivity of the photometric band. If 
we take into account luminosity evolution of the form 



L exp(Qt L /r) 



(13) 



where r is the Hubble time and tr is the lookback time 



t L =r 



f 

J (1 + 2' 



dz' 



{i + z')y/viM{i + z'f + n A 



(14) 



z™ needs to be recalculated accordingly. 

Most of the previous estimates in the literature have de- 
rived the K-correction term by assuming a power-law spec- 
trum of the form S v cx v~ 2 (although in some papers, people 
have tried other methods such as a grey-body model con- 
strained by the 60-100 fim flux ratio). In light of Fig.[T] this 
assumption may not seem very problematic over a narrow 
redshift range except for AGN dust tori. Fig. [4] compares 
the K-corrected 60 fim luminosity derived from the best-fit 
infrared template^ for each galaxy with that derived from a 
S u oc v~~ 2 spectrum over the redshift range 0.003 ^ z ^ 4.0, 
without evolution correction. As expected, the difference be- 
tween the two luminosity estimates is small except for AGNs 
at relatively high redshifts. 

3.2 Evolution 

A similar maximum likelihood approach can be used to de- 
termine the evolution of the LF (S90; Fisher et al. 1992). 
For pure density evolution parametrised in Eq. [6] the evo- 
lutionary rate P can be found by maximising the likelihood 
of observing galaxies at positions rx,r%, ... respectively, that 
is 



n P ( 2l iLo=n 



(1 + ZifdVi 



f z \ (1 + z) p {dV/dz)dz 



(15) 



where z m i n is the lower redshift cutoff. 

The above method, known as the constant-density 
method, is strongly affected by density inhomogeneities in 
the analysed sample. Therefore, a lower redshift limit is usu- 
ally used to reduce this effect which is strongest in the local 
Universe. Table [1] is a compilation of our estimates of P in 
various redshift bins in an Einstein-de Sitter (to facilitate 
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Figure 5. The 1/Vmax 60 fim LF derived for the full sample in 
3 redshift bins without evolution correction. To increase the con- 
trast between LFs in different redshift bins, we have normalised 
the LF by Saunders et al.'s estimate of the present-day LF, i.e. 
Eq. 4. 



comparison with others) or ACDM model and previous es- 
timates. Using the constant-density method, S90 detected a 
strong density evolution P — 6.7 ± 2.3 (Eq. 6). It was con- 
firmed from source number counts in the deep IRAS Faint 
Source Survey (FSS) data base at |6| > 50° (Lonsdale et al. 
1990). However, Fisher et al. (1992) found no evidence for 
evolution in their sample of 5,297 IRAS galaxies flux-limited 
at 1.2 Jy at 60 fim. Oliver et al. (1995) found P « 5.6 in a 
spectroscopic sample of around 2,000 galaxies with Feo ^ 0.2 
Jy. Bertin, Dennefeld & Moshir (1997) also found evidence 
to support strong evolution in a deep subsample of the IRAS 
FSS at a flux limit of Fqo ~ 0.11 Jy, using several methods 
such as background fluctuations (i.e. confusion noise), num- 
ber counts and optical identification rate. SW98 developed a 
minimum variance estimator which is less affected by density 
inhomogeneities and re-analysed the 1.2 Jy samplfl They 
claimed a milder evolution P ~ 4.3 in the IRAS galaxies. 
The IRAS 1 Jy survey of ULIRGs (ultraluminous infrared 
galaxies with L; r 10 12 Lq) showed a strong evolution (Kim 
& Sanders 1998). Takeuchi, Yoshikawa & Ishii (2003) applied 
the constant-density method to the IRAS PSCz and found 
P w 3.4, consistent with SW98. 



3.3 Result and discussion 

Although evolution seems to be detected in most cases, 
there is disagreement in the literature about the magni- 
tude of evolution. The value of P does not vary much 
between an Einstein-de Sitter and ACDM model over the 
limited redshift range examined here. However, it depends 
quite strongly on the redshift boundaries of the sample. Our 



1 For galaxies which are detected at 60 /im only (around 13% of 
the whole sample), we have assumed a power- law spectrum of the 
form Su oc v~ 2 . 



2 The sample used in Fisher et al. (1992) was an early version. 
The final version was used in the re-analysis by SW98. 
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Figure 6. The X/V m3x 60 fim LFs derived for each infrared subpopulation in 3 redshift bins (open circles: 0.02 < z < 0.04; filled circles: 
0.04 < z < 0.10; filled triangles: 0.1 < z < 0.2) without evolution correction. Similar to Fig. 5, we have normalised the LF by Saunders et 
al.'s estimate of the present-day LF (Eq. 4) to increase the contrast between LFs in different bins. Arp 220- type galaxies exhibit strong 
negative evolution. For galaxies with significant AGN contributions, there seems to be some indications of strong positive evolution 
although the error bars are too big for it to be conclusive. Cirrus- and M82-type galaxies show similar weak positive evolutionary trend. 
Because the redshift range is quite narrow, the effect of evolution is not obvious. In Fig. 9, we show that LF derived from maximum 
likelihood method assuming a luminosity evolution of Q = 1.7 (Eq. 13) describes both cirrus and M82 type reasonably well, indicating 
that the two types have similar evolutionary trends. 



adopted value is P — 3.4 derived from the redshift interval 
[0.02, 0.1] in which the sample is expected to be relatively 
complete and less affected by the peculiar velocity field. The 
pure density evolution can be replaced by pure luminosity 
evolution with Q = 1.7 (Eq. 13). In Fig. E] the 1/V max LFs 
of the full sample are shown in three redshift bins within 
redshift z = 0.2. We have corrected for the fact that both 
the first (faintest) and last (brightest) luminosity bins are 
usually biased low. Evolution is clearly seer0 and the mag- 

3 If there is positive evolution, higher redshift points appear 
above the lower redshift points in the luminosity range where 
they overlap. If there is negative evolution, the reverse is true. 
E.g. for density evolution, positive evolution means that number 
density decreases with time (or increases with redshift). 



nitude of evolution is consistent with that derived from the 
maximum likelihood method. 

Galaxy evolution as a function of infrared template type 
is also investigated. Fig. [6]shows the LFs derived for each in- 
frared galaxy population in three bins in the redshift range 
[0.02, 0.2]. Apart from Arp 220-type starbursts, all other 
types seem to experience different degrees of positive evo- 
lutioifl As a result, the dependence of P on the redshift 

4 In Fig. 9, it can be clearly seen that a pure luminosity evolution 
of magnitude Q=1.7 (derived from the whole sample using maxi- 
mum likelihood method) fits the evolution in both cirrus and M82 
subpopulations quite well. AGN dust tori show stronger evolution 
than that of cirrus and M82, while Arp 220 - type galaxies have 
an opposite trend of evolution. 
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Table 1. Estimates of the evolutionary rate P. The last six rows are estimates derived in this paper. 



Sample 

QCD survey (S90) 

The IRAS Faint Source Survey (Lonsdale et al. 1990) 
1.2 Jy survey (Fisher et al. 1992) 
A deep IRAS redshift survey (Oliver et al. 1995) 
A subsample of the IRAS FSS (Bertin et al. 1997) 
1.2 Jy survey (SW98) 

IRAS 1 Jy survey of ULIGs (Kim & Sanders 1998) 
IRAS PSCz (Takeuchi et al. 2003) 

IIFSCz (this paper), 0.01 < z ^ 0.1 (Einstein-de Sitter) 

, 0.02 ^ z < 0.1 (Einstein-de Sitter) 

, 0.02 ^ z < 0.2 (Einstein-de Sitter) 

, 0.01 ^ z ^ 0.1 (ACDM) 

, 0.02 ^ z ^ 0.1 (ACDM) 

, 0.02 ^ z ^ 0.2 (ACDM) 



Method 



6.7 ± 2.3 Constant density method 

~ 7 Source number counts 

2 ± 3 Constant density method 

5.6 ± 2.3 (V/Vmax), number counts 

6.0 ± 1.2 Background fluctuation analysis, number counts 

4.3 ± 2.4 Minimal variance method 

7.6 ± 3.2 Constant density method 

3.4 ± 0.7 Constant density method 

4.7 i 0.9 Constant density method 

3.2 i 0.9 Constant density method 

2.3 i 0.6 Constant density method 

4.3 i 0.9 Constant density method 

3.4 i 0.9 Constant density method 

2.1 i 0.6 Constant density method 




0.20 




0.20 



Figure 7. Left: The logarithmic slope of the selection function versus redshift. The curve is the best-fit analytic SF as parametrised in 
Equation 16. Right: The non-parametric SF (solid line) versus the best-fit analytic SF (dashed line). The normalisation is derived from 
the number density of galaxies with Fgo 0.36 Jy per unit solid angle. 



range is perhaps due to the varying fraction of each infrared 
subpopulation which undergoes different evolutionary trend. 
The apparent negative evolution of the A220-type starbursts 
at low redshift is interesting and perhaps indicates the era of 
major merger events has ended. However, we do not exclude 
the possibility that the very different evolutionary trend seen 
in A220-type starbursts (a very small fraction of the whole 
infrared population) might be caused by incorrect classifica- 
tion. Presumably A220-type starbursts show strong positive 
evolution at higher redshift as a result of higher merger rate. 
Interestingly, this negative evolution at low redshift is pre- 
dicted in the source count model of Rowan- Robinson (2009) . 

The left panel of Fig. [7] displays the local logarithmic 
slope vtih and its errors out to redshift z = 0.2, corrected 
for density evolution P = 3.4. The curve is obtained by a 
minimum \ 2 fitting with the SF parametrised as, 



S(z) 



+ (JL)7]/3/7' 



(16) 



where the best-fit parameters are a = 0.999 ± 0.076, /3 = 
3.587 ± 0.169, 7 = 2.001 ± 0.165, z* = 0.0339 ± 0.0014. 
The right panel of Fig. [7] compares the non-parametric SF 
with the analytic SF. The normalisation is derived from the 
number density of galaxies with Feo zZ 0.36 Jy per unit solid 
angle (~ 2567 sr _1 in the redshift range [0.003, 0.2] in our 
sample) 

1 



N = 



47T 



dz 



(17) 



From Eq. [JJand Eq.[5j the present-day piecewise LF can be 
calculated as 



ML) 



Q (^max) 



Xk 



Sk 



3(-2 m ax)i m i n (2 m ax) 



(18) 



for L min (a;fe-i) < L < L mla (xk)- 

The left panel of Fig. [8] shows LFs derived from vari- 
ous estimators in the redshift range 0.003 ^ z ^ 0.2, after 
correction for a pure luminosity evolution of Q = 1.7. In 
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Figure 9. Luminosity functions at 60 fim derived from the 1/V ma x method (symbols) and the STY method (curves) for cirrus, M82, 
Arp 220 and AGN dust torus type galaxies in the redshift range 0.02 $C z $C 0.1, corrected for a pure luminosity evolution of Q = 1.7. At 
logioC^eo) > 10.3, M82-type starbursts start to dominate in number until about log 10 (I/6o) ~ 11.8. We can also see that the assumed 
luminosity evolution Q = 1.7 describes the LF of cirrus- and M82-type galaxies reasonably well. For AGN dust tori, Q = 1.7 clearly 
underestimates the true magnitude of evolution, while for Arp 220-type starbursts, the degree of evolution is overestimated. 
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determining the LF using the l/Vmux estimator, we weight 
each galaxy by the redshift completeness as a function of 
fro. The errors are calculated from 



E 



weight, 
~V~ 2 ■ 



(19) 



The faint end of the 1/Knax LF is problematic as this 
method is sensitive to density inhomogeneities which are 
more pronounced in the very local Universe. The non- 
parametric LF (derived from the non-parametric SF) shows 
a good agreement with the estimate from S90 and our esti- 
mate using the STY method. Our best-fit parameters in the 
LF parametrised in Eq. [4] are, 



a = 1.29, L* 



W 86r h- 2 L @ , 



0.72. 



(20) 



Therefore, our best-estimate LF using the STY method has 
a steeper faint end than that of S90 and a higher charac- 
teristic luminosity L„. The normalisation is obtained via a 
minimum \ 2 fitting to the 1/Knax LF (squares in Fig. (3), 



X 



E 



(21) 



where (p° bB is the measured 1 /V^nax LF at luminosity Li and 
Oi is the rms error. The normalisation is thus 



C 



,„°bs,„STY /„2 



(22) 



The right panel of Fig. [8] shows the 1/Vmax LF in the 
redshift range 0.2 ^ z ^ 4.0, corrected for a pure luminos- 
ity evolution of Q — 1.7. Again, it is compared with the 
present-day LF of S90. The LF begins to flatten at around 
Lqo > 1O 13 '°L0, which might be due to significant evo- 
lution or a new infrared galaxy population arising at red- 
shift z > 0.2. We have subsequently obtained 99 hours on 
ESO VLT to do spectroscopy for these ultraluminous in- 
frared galaxies (ULIRGs) and hyperluminous infrared galax- 
ies (HLIRGs). The nature of these objects will be discussed 
in future papers. 

The present-day LF for each infrared subpopulation is 
plotted in Fig. [9] assuming a common evolution for all tem- 
plate types. Although Fig.[6]suggests different rates of evolu- 
tion for different spectral types, the samples of the Arp 220- 
type starbursts and AGN dust tori were not large enough to 
determine their rate of evolution reliably, we can see that at 
the low luminosity end, cirrus-type galaxies dominate while 
M82 starbursts dominate at the high luminosity end (at least 
up to 10 12 L Q ). The subpopulation of AGN and Arp 220- 
type starbursts exhibit large deviations from the evolution 
derived for the full sample, reinforcing the conclusion drawn 
from Fig. [6] 



3.4 Monte Carlo analysis of photometric redshift 
errors 

At the flux limit of 0.36 Jy at 60 /im, less than 15% of 
the galaxies in the IIFSCz have photometric redshifts ob- 
tained through the empirical training set (for galaxies with 
near-infrared or radio counterparts) or the template-fitting 
method (for galaxies with optical identifications). Assuming 



the distribution of photometric redshift is a Gaussiarjf], we 
can use Monte Carlo analysis to investigate the uncertainty 
in the LF induced by photometric redshift error. 

We have created 100 realisations of the flux-limited sam- 
ple in the redshift range 0.003 ^ z ^ 0.2 with photomet- 
ric redshifts replaced by random numbers generated from 
a Gaussian distribution {a = 0.02) centred at the original 
photometric redshift estimate, we then derive the 1/Knax 
LF for each realisation. The resulting rms scatter is at most 
7% (mean scatter 4%). Compared to the Poisson error (45% 
at most; mean scatter 10%), we will neglect the uncertainty 
induced by the inaccuracy of the photometric redshifts. 



4 GALAXY BIAS 

4.1 Two-point statistics 

The spatial two-point correlation function £ (r) is one of the 
most frequently used tools in studying galaxy clustering. 
It is defined as the probability of finding a galaxy pair at 
a given separation, in excess of that in a uniform random 
Poisson distribution. In practise, there are several estimators 
available: 



The Blanchard & Alimi (1988; hereafter BA88) estima- 



tor 



f(r) = S^-l 



(23) 



DR(r) n D 

• The Landy & Szalay (1993; hereafter LS93) estimator 



1 



RR(r) 



DD(r)(—f - 2DR(r)(—) + RRM 
no no 



(24) 



• The Hamilton (1993; hereafter H93) estimator, 
DD(r)RR(r) ^ 



£(r) = 



DR 2 (r) 



no and ur are the mean densities of the galaxy and ran- 
dom catalogues respectively. DD(r), DR(r) and RR(r) are 
numbers of weighted galaxy-galaxy pairs, galaxy-random 
pairs and random-random pairs at separation r respectively. 
For volume-limited samples, the weight function applied 
to each galaxy is simply u>i = 1. For flux- limited sam- 
ples, it is common to use the minimum variance weighting 
1/(1 + 4nn(zi) Js(r)), where n(zi) is the density at Zi and 
J'i{ r ) ~ Jo £,{r')r' 2 dr' . The main difference between these 
estimators lies in the sensitivity to the uncertainty of the 
mean density estimation. We found almost identical results 
from these three methods and will only show correlation 
functions obtained from the Landy & Szalay estimator in 
this paper. 

Identical selection effects (both angular and radial) need 
to be reproduced in creating random catalogues without 
structures. We use the luminosity functions of cirrus-type 
and M82-type galaxies to create random samples in calcu- 
lating the clustering of these two subpopulations. As men- 
tioned in Section 2, the redshift completeness of the IIFSCz 
varies across the sky. Therefore, we first apply an angular 
mask which takes into account regions missed by IRAS to 

5 It is a strong (and perhaps not very accurate) assumption that 
photometric redshift errors follow a Gaussian distribution. 
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the generated random catalogue. Then, we modulate the 
number of random galaxies in each grid of the sky accord- 
ing to the redshift completeness so that regions with higher 
redshift completeness will have more random galaxies. 

The well known redshift-space distortions consist of 
small-scale elongation (fingers-of-God effect) and large-scale 
compression of structures along the line-of-sight. The former 
is due to random motions in virialized systems while the lat- 
ter is caused by coherent motions of galaxies towards over- 
dense regions. To derive the real-space correlation function 
£(rV, we calculate the projected two-point correlation func- 
tiorQ (Davis & Peebles 1983) 



p oc 

S((t)=2/ Ar£(«r,7r), 
Jo 



(26) 



where £(cr, 7r) is the two-dimensional correlation function, n 
and a are separations along and perpendicular to the line-of- 
sight respectively. Thus the bias caused by peculiar motions 
along the line-of-sight is integrated out. If the real-space 
two-point correlation function can be well described by a 
power-law £(r) = (r/ro) 1 , then 



(27) 



Here B is the beta function which satisfies 

b[ —. - ) =: - . (28) 



2 ' 2 



r(i) 



To quantify the dependence of galaxy clustering on 
physical parameters such as luminosity and spectral type, 
we define the relative bias between two classes of galaxies as 
a function of scale, 



6i 
62 



(r) 



ft(r) 



(29) 



In this paper, we measure the relative bias of the subsample 
in the form of 



& rc i = [H(cr)/Hfid(o-)] 



1/2 



(30) 



where the assumed fiducial projected correlation function 
corresponds to a real-space correlation function f(r) = 
(r/3.5)" 1 ' 8 . The same approach was used in Zehavi et al. 
(2005). 



4.2 Error analysis 

Due to the correlated nature of the correlation function at 
different separations, Poisson error may seriously underesti- 
mate the true uncertainties in the correlation function. To 
determine the covariance of our data, we use the bootstrap 
resampling technique (Ling, Frenk & Barrow 1986) to create 
a series of data sets. 



A( Si ) = 



cr(sj) 



(31) 



6 The projected correlation function is sometimes denoted as 
w p (r p ), where r p is the separation perpendicular to the line-of- 
sight. 



where f(si) is the mean value of the correlation function at 
separation Si from the different bootstrap samples and <r(sj) 
is the rms scatter. The covariance error matrix is 



Ca = (A( Si )A( Sj )}. 



(32) 



The maximum-likelihood estimate of the parameters in the 
model can be calculated by minimising 

X 2 = (A rcs ) T C" 1 A rcs , (33) 

where A rcs = (£° b — cd _ £ modcl )/cr. However, the boot- 
strap resampling method does not account for uncertainties 
caused by fluctuations outside the sampled volume. In all 
correlation functions presented below, we have used the full 
covariance matrix to derive the best-fitting power law. 



4.3 The luminosity dependence of clustering 

The infrared emission of galaxies comes from dust emission 
associated with regions of active star formation as well as 
the general interstellar radiation field. Luminous infrared 
galaxies are often starbursts or dust-enshrouded AGN. As 
a result, the FIR luminosity is expected to be an excel- 
lent SFR tracer for strong compact dusty starbursts, but 
is less reliable for the disks of normal galaxies where radi- 
ation from old stars contribute significantly to the heating 
of dust (Helou 1986; Kennicutt 1988; Rowan-Robinson et 
al. 1997). In this paper, we use the 60 micron luminosity to 
derive SFR (Rowan- Robinson et al. 1997), 



(34) 



where e = 2/3 is the adopted fraction of UV light absorbed 
by dust. 

Previous studies using smaller samples were unable to 
detect any luminosity-dependence of the clustering ampli- 
tude, possibly due to the contrasting effects of the dark 
matter halo mass and specific star formation rate (SSFR) 
on determining the far-infrared luminositjQ Based upon the 
luminosity - colour correlation and the colour dependence of 
the correlation function, Hawkins et al. (2001) proposed a 
luminosity-dependence of the form 



log 10 (s ) = -0.028 log 10 (L 60 ) + 0.95, 



(35) 



where so is the correlation length in the redshift space. So, 
unlike optical galaxies, they found a weaker correlation for 
infrared galaxies of higher luminosity but their result is of 
low statistical significance. However, there might be a prob- 
lem with their estimate because they used the redshift-space 
correlation function. The effect of redshift distortion at the 
scale we are interested in tends to boost the clustering signal 
and it is worse for nearby samples where peculiar velocities 



7 The SSFR is lower in more massive galaxies at both the lo- 
cal Universe z ~ and the distant Universe z ~ 1 (Brinch- 
mann et al. 2004; Elbaz et al. 2007; Noeske et al. 2007; Zheng 
et al. 2007; Damen et al. 2009). Oliver et al. (in preparation) 
use the Spitzer Wide-area Infrared Extragalactic Legacy Sur- 
vey (SWIRE) to study the SSFR as a function of stellar mass 
parametrised as SSFR oc , where M, is stellar mass. They 
find that the SSFR shows a less steep decline for late-type galax- 
ies (/3 ~ —0.15) than early type galaxies (J3 ~ —0.46). 
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Table 2. Nested volume-limited samples of all galaxies (from VI to V7), M82-type starbursts (from MV1 to MV6) and cirrus-type 
quiescent galaxies (from CV1 to CV5). The columns are sample name, redshift range, 60 micron luminosity threshold, mean 60 fim 
luminosity, number of galaxies, correlation length rrj, power-law index 7, mean bolometric infrared luminosity derived from the best-fit 
infrared template and SFR. 



Sample z range 



Min(log(L 60 )) <log(L 60 )> 



r (/i" 1 Mpc) 7 



(log(L IR )) (SFR) (M Q /yr) 



VI 
V2 
V3 
VI 
V5 
V6 
V7 



[0.003, 0.01] 
[0.003, 0.02] 
[0.003, 0.03] 
[0.003, 0.04] 
[0.003, 0.05] 
[0.003, 0.06] 
[0.003, 0.07] 



9.00 

9.60 

9.96 

10.22 

10.42 

10.58 

10.72 



9.48 ± 0.38 
9.94 ±0.28 
10.24 ±0.24 
10.46 ± 0.22 
10.64 ±0.21 
10.80 ±0.20 
10.92 ±0.19 



1010 
2352 
3081 
3283 
2916 
2617 
2349 



3.75 ± 0.46 
3.30 ±0.21 
2.85 ±0.24 
3.42 ±0.35 
3.50 ±0.43 
3.05 ±0.45 
3.03 ± 1.08 



2.55 ± 0.24 
1.67 ±0.06 
1.69 ± 0.07 
1.74 ±0.09 
1.66 ± 0.10 
1.60 ±0.16 
1.83 ±0.37 



9.91 ±0.37 
10.33 ± 0.29 
10.55 ± 0.26 
10.79 ± 0.29 
11.00 ±0.25 
11.10 ±0.24 
11.26 ±0.23 



1.00 

2.87 

5.71 

9.56 

14.50 

20.76 

27.68 



Sample z range 



Min(log(L 60 )) (log(L 60 )> N 3 



r (h 72 Mpc) 7 



(log(L 7Jl )> (SFR) (M Q /yr) 



MV1 
MV2 
MV3 
MV4 
MV5 
MV6 



[0.003, 0.01] 
[0.003, 0.02] 
[0.003, 0.03] 
[0.003, 0.04] 
[0.003, 0.05] 
[0.003, 0.06] 



9.00 

9.60 

9.96 

10.22 

10.42 

10.58 



9.61 ±0.44 
10.01 ± 0.31 
10.29 ± 0.26 
10.49 ± 0.23 
10.66 ± 0.22 
10.81 ± 0.20 



343 

1050 

1650 

1984 

1990 

1880 



3.39 ± 1.20 
2.75 ±0.39 
2.60 ±0.83 
3.54 ±0.46 
3.38 ±0.87 
3.99 ± 1.08 



1.89 ± 0.42 

1.63 ± 0.14 

1.64 ±0.24 
1.89 ± 0.16 
1.64 ±0.20 
2.01 ±0.36 



9.97 ±0.44 
10.36 ± 0.30 
10.57 ±0.26 
10.80 ±0.30 
11.01 ± 0.25 
11.10 ±0.23 



1.36 

3.40 

6.42 

10.30 

15.21 

21.13 



Sample z range 



Min(log(L 60 )) (log(L 60 )> N 9 



r (h T2 Mpc) 7 



<log(Li fl )> (SFR) (M Q /yr) 



CV1 
CV2 
CV3 
CV4 
CV5 



[0.003, 0.01] 
[0.003, 0.02] 
[0.003, 0.03] 
[0.003, 0.04] 
[0.003, 0.05] 



9.00 

9.60 

9.96 

10.22 

10.42 



9.42 ± 0.32 
9.88 ± 0.23 
10.17 ±0.18 
10.39 ± 0.15 
10.56 ±0.12 



624 

1158 

1157 

972 

580 



3.56 ±0.43 
4.27 ±0.20 
2.41 ± 1.01 
4.31 ± 1.12 
3.47 ±2.16 



2.41 ± 0.50 
1.87 ±0.06 
1.51 ± 0.22 
1.84 ±0.26 
1.69 ± 1.02 



9.89 ±0.32 
10.34 ±0.22 
10.56 ± 0.19 
10.80 ±0.25 
11.01 ±0.21 



0.86 
2.52 
4.93 
8.12 
11.99 



along the line-of-sight are not negligible. Thus, the correla- 
tion length derived from redshift-space correlation function 
will be higher than the true value for nearby samples. It is 
also important to bear in mind that Eq. 35 was proposed 
to describe the luminosity-dependence for the full sample 
including both cirrus and M82 type. Therefore, it is not re- 
ally relevant to the question we are trying to address in this 
paper, i.e. the relation between SFR and environment for 
star-forming galaxies. 

To investigate the luminosity dependence of galaxy clus- 
tering and at what range it occurs, the IIFSCz is divided 
into nested volume-limited subsamples. In Table [2J we have 
listed redshift range, 60 /im luminosity threshold, mean 60 
(jm luminosity, number of galaxies in the sample, correlation 
length ro in the real-space, power-law index 7, the bolo- 
metric infrared luminosity derived from the best-fit infrared 
template and SFR derived from Eq. 34. We also extract 
volume- limited samples of cirrus- type galaxies and M82-type 
starbursts. In Fig. 1101 the relative bias defined by Eq. 30 with 
(j = 3 Mpc is plotted against the mean 60 micron luminos- 
ity Lqo for each subsample in Table The error bars are 
the rms scatter derived from 25 bootstrap realisations. The 
empty stars and the filled stars represent the relative bias 
for cirrus and M82 subsamples respectively. The open cir- 
cles show the relative bias of the full sample, which closely 
follow the relative bias of cirrus-type galaxies at low lumi- 
nosities log 10 (L6o) < 10.3 and then M82-type starbursts at 
higher luminosities. This behaviour of the whole sample is 
expected as we have already seen from Fig. 3 and Fig. 9 
that cirrus-type galaxies is the dominant subpopulation at 
low redshift / luminosity and then M82-type starbursts take 
over at around log 10 (Z/6o) ~ 10.3. 



For cirrus-type quiescent galaxies, dust emission caused 
by old stellar population contributes significantly to the in- 
frared emission and thus the derived SFR are subject to 
large uncertainties. To study the relation between current 
star formation activity and environment, it is more relevant 
to look at the clustering of M82-type starbursts as a function 
of infrared luminosity. In Fig. 1101 we can see that the cor- 
relation strength correlates tightly with SFR for M82-type 
starbursts. The thick solid line is the best-fitting straight 
line of the SFR - clustering strength correlation. Possible 
cause for this tight correlation will be discussed in Section 
4.5 where our result is compared with clustering of star- 
forming galaxies in the distance Universe (at 2 ~ 1). 

4.4 The infrared template type dependence of 
clustering 

The top panel of Fig. 1 in Wang & Rowan- Robinson (2008) 
shows the ratio of -F100 and Fqo as a function of spectro- 
scopic redshift compared with predictions based on four in- 
frared templates. At redshift 2 < 1.2, cirrus-type quiescent 
galaxies generally have cooler FIR colours than M82 star- 
bursts. In Hawkins et al. (2001), cold galaxies are defined by 
-Fioo/^60 > 2.3 which are essentially cirrus-type galaxies in 
our sample. The colour-luminosity correlation means galax- 
ies of different colour will preferentially be found in different 
volumes. To reduce the effect of cosmic variance and lumi- 
nosity dependence of clustering (see the previous section), 
we need to compare subsamples of cirrus-type galaxies and 
M82 starbursts, corresponding to cooler and hotter galaxies 
respectively, in the same volume (redshift range) and same 
luminosity bin. 
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Figure 10. Relative bias for volume-limited samples listed in Table 2 as a function of 60 micron luminosity. The relative bias is calculated 
using Eg. 1301 and a = 3 Mpc. The error bars are derived from 25 bootstrap realisations of each subsample. At log 10 (L6o) < 10.3, cirrus- 
type galaxies dominate in number (see Fig. 9), so galaxy clustering of the whole sample (open circles) follows that of cirrus-type (empty 
stars). At higher luminosities log 10 (L6o) > 10.3, M82-type starbursts start to dominate and the clustering of the full sample is almost 
the same as that of M82-type starbursts. The thick solid line is the best-fit straight line of the luminosity-dependence of clustering for 
M82-type starbursts which clearly exhibits a correlation between SFR and clustering strength. 



Table 3. Flux-limited subsamples of M82-type starbursts (from MFl to MF3) and cirrus-type galaxies (from CFl to CF3). The columns 
are sample name, redshift range, mean 60 fim luminosity, number of galaxies, correlation length and the power-law slope. 



Sample 


z range 


(log(L 60 )) 




r (h 7 J Mpc) 


7 


MFl 
MF2 
MF3 


[0.003,0.015] 

[0.015,0.03] 

[0.03,0.05] 


9.52 ±0.55 
10.11 ±0.33 
10.49 ±0.27 


1027 
2205 
2777 


3.52 ±0.31 
2.65 ±0.50 
3.96 ± 0.30 


1.73 ±0.11 
1.62 ±0.13 
1.76 ±0.10 


Sample 


z range 


(log(i60)) 


Kg 


r {h~^ Mpc) 


7 


CFl 
CF2 
CF3 


[0.003,0.015] 

[0.015,0.03] 

[0.03,0.05] 


9.34 ± 0.44 
9.95 ±0.26 
10.35 ±0.19 


1710 
2225 
1440 


4.05 ±0.20 
3.15 ±0.23 
4.93 ±0.31 


1.78 ±0.05 
1.74 ±0.07 
2.02 ±0.17 



Already in Fig. 10, we can see that cirrus- type galaxies 
have larger relative bias than M82-type starburst in all lu- 
minosity bins. Fig. [TTJcompares the projected two-point cor- 
relations functions of the two infrared subpopulations. The 
top row shows E(<t)/<7 for selected volume-limited samples of 
the two types listed in Table 2. The selection is based on the 
number of galaxies to ensure that we are not comparing very 
noisy correlation functions. Similarly, the bottom row shows 
the H(a) jo for selected flux-limited samples of the two types 



listed in Table 3. Difference in clustering strength between 
cirrus and M82 is clearly present in the sense that cirrus- type 
galaxies are more strongly clustered than M82 starbursts. 
We can also derive the correlation lengths with fixed power- 
law index y = 1.8 for subsamples of cirrus and M82-type 
galaxies listed in Table 2 and the relative bias is calculated to 
be b c irrus/&M82 = 1-25 ± 0.07 on scales 2 ^ r ^ 13 Mpc. 
At scales larger than ~10 Mpc, there is some indication that 
the relative bias between the two infrared types decreases. 
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Figure 11. Top row: The projected two-point correlation function for selected volume-limited samples of cirrus-type (filled circles) and 
M82-type (open circles) galaxies listed in Table 2. The dashed (dotted) line is the best-fitting power law to the cirrus (M82) correlation 
function on scales 2 ^ a ^ 13 Mpc, using the full covariance matrix. The error bars represent the rms spread derived from 25 bootstrap 
realisations. Bottom row: The projected two-point correlation function for selected flux-limited samples of cirrus-type and M82-type 
galaxies listed in Tabic 3. 



Our result is consistent with the conclusions of Hawkins et 
al. (2001) that cooler galaxies cluster more strongly than 
hotter galaxies. 

The clustering dependence on infrared template type is 
perhaps analogous to the colour (or optical spectral type) 
dependence found in optical galaxies (Madgwick et al. 2003; 
Zehavi et al. 2002, 2005; Wang et al. 2007; Cresswell & Perci- 
val 2009; Ma et al. 2009). It suggests that the bulk of infrared 
luminosity for the two types, cirrus and M82, arises from 
different sources. The FIR emission of cirrus-type galaxies 
could be mainly due to dust heated by the general interstel- 
lar radiation field where old stars dominate. Given the local 
SFR - density relation, star-forming galaxies are preferen- 
tially found in low-density regions and so it is expected that 
the clustering strength of M82-like starbursts is weaker than 
that of cirrus-type quiescent galaxies. 

4.5 Comparison with the distant Universe (z ~ 1) 

At z ~ 1 , Cooper et al. (2006) found luminous blue galaxies 
in the DEEP2 Galaxy Redshift Survey are found in environ- 



ments with greater overdensity. There are also evidence that 
clustering strength of star-forming galaxies correlates with 
SFR in the distant Universe. Gilli et al. (2007) used star- 
forming galaxies detected at 24 micron by Spitzer/MIPS 
in the Great Observatories Origins Deep Survey (GOODS) 
fields to measure the spatial clustering. They found galaxies 
of higher infrared luminosity have a larger correlation length 
than less luminous ones. It indicates that more intensive star 
formation is hosted by more massive dark matter halos. 

The stronger clustering of star-forming galaxies with 
higher SFR detected in the distant Universe can be ex- 
plained by the tight correlation between SFR and stellar 
mass (Noeske et al. 2007; Elbaz et al. 2007), assuming more 
massive galaxies generally reside in more massive dark mat- 
ter halos. However, because the SFR - stellar mass relation 
is also valid in the local Universe, we should expect a sim- 
ilar correlation between clustering strength and SFR which 
is exactly what we see in the clustering dependence on in- 
frared luminosity for M82-type starbursts (the solid line in 
Fig. 10). 
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5 DISCUSSION AND CONCLUSION 

We have derived the luminosity function and selection func- 
tion of the IIFSCz to learn about the evolution of infrared 
galaxies. To minimize the effect of incompleteness, we chose 
a flux limit of 0.36 Jy at 60 fim. Either a pure density evo- 
lution P — 3.4 or pure luminosity evolution Q — 1.7, is 
detected in the full sample using the maximum likelihood 
method. The magnitude of evolution is at the lower end of 
previous estimates. The value of P shows some dependence 
on the redshift boundaries which might be caused by dif- 
ferent mixtures of infrared subpopulations which have very 
different evolutionary trends. The present-epoch 60 (jm lu- 
minosity function derived from galaxies in the redshift range 
0.003 < z < 0.2 shows good agreement with the luminosity 
function from Saunders et al. (1990). However, the lumi- 
nosity function derived from galaxies in the redshift range 
0.2 < z < 4.0 shows evidence for much stronger evolution 
or possibly an entirely new class of infrared galaxies of high 
luminosity emerging at redshift z > 0.2. Follow-up spec- 
troscopy studies will investigate further the population at 
the high luminosity end. 

The spatial clustering of galaxies as a function of galaxy 
properties such as type and luminosity provides strong con- 
straints on models of galaxy formation and evolution. In 
particular, the clustering of the infrared galaxies contains 
essential information on the relation between star forma- 
tion activity and environment. We use the projected corre- 
lation function H(<r), which is free from redshift-space dis- 
tortion, to study the clustering dependence on infrared tem- 
plate type and infrared luminosity. When comparing clus- 
tering strength between the two infrared subpopulations, 
we found that cirrus-type galaxies cluster more strongly 
than M82 starbursts in all luminosity bins until the cirrus- 
population becomes absent at high luminosity /redshift. The 
average relative bias between the two infrared populations is 
&cirrus/&M82 = 1.25±0.07 on scales 2 r sC 13 Mpc. This 
is consistent with the SFR - density relation which states 
that star-forming galaxies preferentially reside in low den- 
sity environments. The dependence of clustering strength on 
SFR (using infrared luminosity as a proxy) is investigated 
using M82-type starbursts. We detect a correlation between 
SFR and clustering strength, which follows the basic trend 
found in the distant Universe (at z ~ 1). It can probably be 
explained by the tight correlation between SFR and stellar 
mass which holds at the local Universe as well as z ~ 1 and 
possibly z ~ 2. However, the trend that the SFR increases 
with clustering strength must becomes invalid at some point, 
as we do not find galaxies with the highest SFR in environ- 
ments of highest density (such as centres of rich clusters) 
in the local Universe. The two seemingly inconsistent state- 
ments, the SFR - density relation and the correlation be- 
tween SFR and clustering strength, can be reconciled if the 
latter operates mainly in the low density regime and be- 
comes invalid once the density exceeds a certain threshold. 

Direct measurement of the infrared emission of galaxies 
is crucial to measuring SFR of distant galaxies. We hope to 
extend our study of the relation between SFR and environ- 
ment to higher redshift using surveys such as SWIRE and 
COSMOS. Observations to be made by the Herschel satellite 
will undoubtedly provide more insight. 
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